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Abstract 

Jamming is a phenomenon occurring in systems as diverse as traffic, colloidal suspensions and 
granular materials. A theory on the reversible elastic deformation of jammed states is presented. 
First, an explicit granular stress-strain relation is derived that captures many relevant features of 
sand, including especially the Coulomb yield surface and a third-order jamming transition. Then 
this approach is generalized, and employed to consider jammed magneto- and electro-rheological 
fluids, again producing results that compare well to experiments and simulations. 
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We start our study of jamming [1] in granular systems, by deriving an appropriate stress- 
strain relation from a simple, postulated elastic energy. It accounts for the reversible elastic 
deformation of granular systems, up to the point of vield, and reproduces many relevant 
results from granular physics and soil mechanics 0, y], including the compliance tensor, 
Rankine states, and shear dilatancy. Moreover, the elastic energy is convex only below the 
Coulomb yield condition and becomes unstable there. As a result, the system escapes from 
the strained state and looses shape-rigidity, providing an explanation why sand unjams. 
Next, the granular elastic energy is shown to be a special case of a more generally valid 
energy expansion, with respect to the shear strain. Realizing that this expansion may serve 
as the starting point to account for other jammed systems, we use it to consider colloidal 
suspensions j^|, specifically magneto- and electro-rheological fluids, which solidify at fields 
strong enough Q]. Again, an energy expression is proposed, from which the magnetic, 
dielectric and elastic behavior is deduced, especially the solid-fluid phase diagram. 

Our basic understanding of sand is due to Coulomb, who noted that its most conspicuous 
property is yield: A pile of dry sand possesses a critical slope that it will not exceed. His 
insightful conclusion is that the quotient of shear stress over pressure must not exceed a 
certain value, \o~ s \/P < fif. Wet sand can sustain a small shear stress ov even at vanishing 
pressure. It satisfies the Mohr-Coulomb condition, |er s | < fijP + <j c , see |7[. 

It is standard praxis in soil mechanics to calculate the stress distribution by taking the 
stress <Jij as some function of the strain Uij. Unfortunately, the calculated stress distribution 
routinely contradicts the Coulomb condition, and yield must be postulated, ex post facto, 
where it is not satisfied. An improvement of this somewhat brute method is given by the 
Rankine states, a s = ±P[if, which should hold close to yield. The ameliorated calculation 
is given by accepting the result of elasticity away from the region of failure, postulating a 
Rankine state close to it, and connecting both smoothly. Clearly, in spite of ingenious ways 
to circumvent it, the basic problem is the lack of a stress-strain relation Uij(o~ki), with which 
a realistic stress distribution can be calculated. 

If we had Uij(o~ki), the incremental relation, 5uij = (duij/daki)Saki = Xiju^Cki, is easily 
derived. The elements of the compliance tensor \jki can also be obtained from experiments, 
in which 5uij, the strain response to a stress change 5o~ij, is measured |8[. Although integrat- 
ing the measured \ijki should in principle lead to Uij(aki), this is a hard, backward operation 
- made more difficult by the typical scatter of data, partly from irreversible plastic deforma- 
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tions. This circumstance has led many to espouse the view that \ijki is history-dependent, 
that an explicit Uij((Tki) (from which to deduce Xijki) does not exist. Different elasto-plastic 
theories, some exceedingly complex, have been constructed to account for \jku including 
both elastic and plastic deformations, though a universally accepted model is missing j^. 

Confining our study to reversible elastic deformations, we derive a stress-strain relation 
to account for the listed granular behavior. We start from the elastic energy 

w = l5°- 5 (B5 2 + Au 2 a ), (1) 

where 5 = — uu is the compression, u 2 = u°jU°j is shear strain squared, (uu denotes the 
trace of the strain and u°j its traceless part. 5, u s = imply the grains are in contact but 
not compressed or sheared.) A, B > are functions of the void ratio e, an independent 
variable. We adopt the same empirical expression for both, A, B ~ (2.17 — e) 2 /(l + e), 
see js . Eq (JJ) is clearly evocative of the Hertz contact: The energy of compressing two 
elastic spheres scales with (Ah) 2 ' 5 , where Ah is the change in height j9(- Writing the energy 
as ^E(Ah) 2 , the effective Young modulus E ~ (Ah) ' 5 vanishes with Ah. The physics for 
the shear modulus is assumed to be similar. 

We postulate Eq (0) to consider its ramifications - noting that it should be possible to 
derive it employing micro-mechanics n]]: Although an intricate task, it is not as difficult as 
calculating the stress or the compliance tensor \jki directly. Remarkably, assuming that 
both moduli vanish with 5 0,5 , we take sand to be arbitrarily pliable, not at all "fragile" 11 1 . 
Differentiating the energy w with respect to 5, u s yields the pressure P and shear cr s , two 
scaler quantities; differentiating it with respect to Uy yields the complete stress tensor cx^, 

P = dw/05 = § B 5 1 - 5 + \A u 2 J5 - 5 , (2) 
a s = dw/du s = A5°- 5 u s . (3) 
(Tij = dw/duij = -P5 i:j + A5°- 5 ul. (4) 

This is the announced static stress-strain relation. The first term in P is well-known and 
considered characteristic of Hertz contacts. The second term, accounting both for shear 
dilatancy and yield, is new. Dilatancy: Holding P constant, 5 decreases (and the volume 
expands) with growing u s . Yield: For given u s , the compressibility (dP /d5)~ l is negative 
if 5 is sufficiently small. This implies lack of local stability, and the system will not re- 
main in the strained state. It is then, without the capability to sustain static shear, in 



a fundamental sense "fluid" - though by no means necessarily Newtonian. In fact, the 
energy looses stability even before dP/d5 turns negative, as the cross convexity condition 
(d 2 w / dd~ 2 )(d 2 w / du 2 ) > (d 2 w/d5 du s ) 2 , or u 2 / S 2 < 5B/2A, also needs to be met. We saw the 
significance of instability in a previous work [l2[ , but did not realize the following remarkable 
point and its consequences: Rewriting the cross convexity condition by replacing 5, u s with 
P, a s leads directly to (the Drucker-Prager variant [7( of) the Coulomb yield condition, 

\a s \/P<^/4A/5B. (5) 

To account for wet sand, the term — P C S (with P c > 0) is added to the energy w. This implies 
a force (typically supplied by the water's surface tension) that compresses the grains even 
without an applied pressure. The additional term does not change the convexity condition, 
only substitutes P + P c for P in Eq (J2J). As a result, Eq (j5J) assumes the Mohr- Coulomb 
form, \a s \ < (P + P c ) ^/4A/5B. 

As any other elasticity theory, the stress-strain relation of Eqs (I2I-]I4J) may be directly 
solved with appropriate boundary conditions to obtain a complete stress distribution. Be- 
cause it includes yield as given by Eq (jSJ), the Rankine states are automatically predicted 
where instability is close. And the compliance tensor \ijki is obtained by simple differenti- 
ation. Writing Suij = XijkiSo'ki as a vector equation, 5a = M5u, with M a 6 x 6 matrix, 
we see yield is signified if an Eigenvalue mi of M vanishes, with the Eigenvector 5u\ indi- 
cating the direction of instability. Explicit calculation shows 5ui\\(dmi/da), implying Sui 
is perpendicular to the yield surface, mi(a) ~ \o~ s | - P^/AA/5B = 0. If there is no plastic 
contribution, this implies flows perpendicular to the yield surface, a circumstance referred 
to as the "associated flow rule" |2(. 

In view of these results, there can be little doubt that Eq (Q) indeed captures the essence 
of granular elasticity. And the remaining question is: To which extent is it also a quantitative 
rendition. To test this, we compare the calculated Xijki to the data gathered recently 
over a wide range of pressure, shear stress and void ratio. (Specifying these three variables, 
the reversible granular response is unique, showing no history-dependence.) Fig. ^ is a 
typical plot, with an overall agreement that further confirms Eq ([T]) . (The expression for 
Xijki is too cumbersome to be displayed here. It will be given in a forthcoming single- 
issue paper containing extensive comparison.) Note the ratio A/B is fixed by the Coulomb 
friction coefficient /i/, so the theory has only one overall scale factor, and no actual adjustable 
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FIG. 1: The Poisson ratios v zx , v xy , the Young moduli E Z ,E X = E y , and the shear moduli G zx , G xy , 
measured [a] with Ham River sand at cr xx /a zz = 0.45 and a void ratio of 0.66, compared to the 
calculated curves assuming B = = 6800Mpa, with E{ = A^, G%j = Vij = —Xajj/Xuu- 

(x, y are horizontal directions, z the vertical one.) Note these coefficients are pairwise equal for 
linear elasticity, but deviate from each other nonlinearly; theory and experiment especially agree 
with respect to the direction of deviations, ie., the fact that v zx > u xy , E z > E x = E y , G zx = G xy . 



parameter. (The most important effect missing in Eq JI} is probably "fabric-anisotropy" |2|.) 

Switching now to a broader context, we proceed to discriminate between the general 
feature of the above theory and those aspects specific to granular elasticity. This should 
give us a better appreciation why Eq ((TJ) is as successful, and also help to apply the same 
approach to other jammed systems. Generally speaking, the energy should be a function of 
at least two variables, u s and /, with / being the one driving the transition, taking place 
at f c . In sand, suspensions, and electro-rheological fluids, / is respectively given by the 
compression 5, concentration, and the electric field. Expanding the energy in u s , 

w = w (f) + \Kul (6) 

the shear modulus K is a function of /, typically K ~ (/ — f c ) a with a > in the solid 
phase (/ > L), and K = in the liquid one (/ < f c ). This dependence is observed in 
suspensions j^, simulations [3| and, with a ~ |, works well for sand. We take it as an 
input. Local stability requires K > and 

wl>[{K'flK-\K"]ul^Kul (7) 



ensuring w is convex in f,u s . Because k ~ a(a + 1) x (f c — f) a ~ 2 is positive, the inequality 
is always violated when u s becomes sufficiently large, rendering instability, and hence the 
unjamming transition, a generic feature. If a < 2, k diverges for / — > f c , and unjamming 
occurs at vanishing values of u s (assuming w' Q ' remains finite). This ensures the validity of 
the expansion of Eq (jHJ) . 

Considering the jamming transition in the shear- free limit u s — > 0, we identify it - by 
analogy to conventional phase transitions - as of n th order, if d l wo/df' 1 is continuous for 
i < n, but not for i = n. With u> ~ S 25 , sand displays a third-order jamming transition. 

Yield at finite shear, as a result of the energetic instability, Eq (J2J), is not an equilibrium 
transition, because the liquid phase moves and dissipates. This may well be compared 
to raising the temperature T in a current-carrying superconductor, such that the metal is 
pushed into its normal state carrying a dissipative, ohmic current. In fact, if one identifies / 
as T, replaces u s with the superfluid velocity v s (and hence a s with the current, j s = p s v s ), 
Eq (JHJ) is valid for superconductors, and superfluid helium, \Ku 2 s — > \p s v 2 s , respectively with 
p s ~ T c — T and p s ~ (T c — T) 2 / 3 |l5| . Macroscopically, jamming and phase transition are 
clearly hard to tell apart, and their conceptual difference must be subtle. 

Next, we consider ER and MR (or electro- and magneto-rheological) fluids, employing 
them as further examples for the above notion of jamming. Although experimental data are 
as yet not confining enough for an unambiguous determination of their energy, plausibility 
may be drawn on to fill the gap. In ER fluids, the dielectric displacement D assumes the role 
of the transition-driving variable /. Writing the shear- free part of the energy as wq = W1+W2, 
we take W\ = ^D 2 /ei, accounting for a linear dielectric relation, and W2 = — |A(Z? — D c ) 2 , 
assuming that linearity prevails after the transition at D c . This is a second-order transition, 
and the electric field E = dw/dD has a kink at D c : We have E = D/e\ for D < D c , and 
E — E c = (D — D c )/e 2 for D > D c , with l/e 2 = 1/ei - A, E c = D c jt x . (A discontinuity in 
E, or a first-order transition, was to our knowledge never reported. Higher order transitions 
are possible, seem even likely, but they are not compatible with a linear constitutive relation 
after the transition.) u>2 is the condensation energy, so A must be positive for solidification 
to take place. (Taking D — > B, E —>■ H yields the analogous formulas for MR fluids.) Given 
w and K = A(D — D c ) a , the energy a; of Eq © is specified. We calculate the dielectric 
relation E = dw/dD\ Ug , elastic relation a s = dw /8u s \d, and rewrite Eq (JJJ) in terms of E, cr s 
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FIG. 2: Elastic and dielectric properties, including the yield point, for electro- and magneto- 
rheological fluids: Shear stress a s versus shear strain u s at fixed electric field E, and E — E c versus 
D — D c (or H — H c versus B — B c ) at fixed cr s , for the exponents a = 2, 1,0.2. Choosing the 
dimension of both curves such that the yield points are at (1,1) render the curves universal - 
removing the dependency on (i) all material parameters other than the exponent a, (ii) E in the 
upper plot, and (iii) a s in the lower one. The relation a s (u s ) for a = 1,2 agree with data from 
experiments and simulations jlgl ]. 

to obtain the yield condition, 



The exponent a = 1, or a yield stress |<7 S | ~ (H — H c ) 3 ^ 2 is observed for most MR-fluids [l 



The same value is also appropriate for a few ER-fluids , though the yield stress is typically 
quadratic j^Q], \cr s \ ~ (E — E c ) 2 , indicating a = 2. An ER-fmid capable of sustaining an 
unusually high shear strength was reported p to display a nearly linear dependence of the 
yield stress, \a s \ ~ E — E c , or a < 1. For a = 0, the shear modulus K is independent of the 
field, and there is no yield at all. This is the reason the square root in Eq (jHJ) diverges for 
a — > 0, and possibly explains the observed high yield stress. 

Finally, the above approach and results are critically appraised. (Granular vocabulary is 
employed for this purpose, though the statements are equally valid for ER and MR fluids.) 
In physics, every microscopic state has a unique energy. The same holds for macroscopic ones 
if we insist on a consistent description. The macroscopic energy always depends on entropy 
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and conserved quantities, such as momentum and mass density. And if the considered system 
can sustain static shear stresses, the strain field Uij must also be included as an independent 
variable, where is to be understood, in soil-mechanical parlance, as the reversible elastic 
portion of the strain field. 

It is a plain fact that sand piles, if left alone under gravity, are stable - in spite of every 
kind of infinitesimal perturbations, which are always present. This demonstrates sand's 
capability to sustain static shear and is the reason for including u^j. Irrespective whether a 
unique displacement field exists, the elastic description employing is robust enough to be 
valid. This is not different from superfluid helium with vortex lines, in which the description 
in terms of the velocity v s = — V0 remains sound, although the phase field is multivalued. 

Given an energy expression w, its derivative dw/duij yields the stress tensor cr^-, and its 
second derivative d 2 w / 'du^du^i the inverse of the compliance tensor In soil mechanics, 
the usual approach consists of postulating the stress dependence of the 18 independent 
components of \ijki directly, while seeking the account for the plastic contribution at the 
same time, referring to the result as constitutive relations [3]. This is quite obviously a much 
harder task than finding the one appropriate scalar expression for the energy w which, even 
if heavy-handedly simplified, preserves a large number of geometric correlation by the mere 
fact that \jki is obtained via a double differentiation. We believe this to be the main reason 
why the calculated Ayjy stood up so well when compared to the extensive data of ^]. 

The expression we proposed in Eq (|T} is indeed the result of weighing simplicity versus 
accuracy while stressing the former, and hence is subject to further scrutiny. As discussed, 
it includes first of all an exp ansion in u s : w = \Ku 2 s assuming K ~ 5 a . Starting from 
w = \{B8 2+b + A8 a u 2 ) in |l2j], we considered the experiments of inclined plane, simple shear 
and triaxial test to arrive at a ~ 0.4, b ~ 0.5 giving the best agreement On the other 
hand is the fact that the Coulomb yield condition, Eq (|5|). remains unchanged as long as 
a = b. And it becomes implicit if a, b deviate from each other - though the numerical 
difference is at first modest. Our tentative choice is a — b — ~. 
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